Multi-fluid cosmology: An illustration of fundamental principles 
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Our current understanding of the Universe depends on the interplay of several distinct "matter" components, 
which interact mainly through gravity, and electromagnetic radiation. The nature of the different components, 
and possible interactions, tends to be based on the notion of coupled perfect fluids (or scalar fields). This 
approach is somewhat naive, especially if one wants to be able to consider issues involving heat flow, dissipative 
mechanisms, or Bose-Einstein condensation of dark matter We argue that a more natural starting point would 
be the multi-purpose variational relativistic multi-fluid system that has so far mainly been applied to neutron 
star astrophysics. As an illustration of the fundamental principles involved, we develop the formalism for 
determining the non-linear cosmological solutions to the Einstein equations for a general relativistic two-fluid 
model for a coupled system of matter (non-zero rest mass) and "radiation" (zero rest mass). The two fluids 
are allowed to interpenetrate and exhibit a relative flow with respect to each other, implying, in general, an 
anisotropic Universe. We use initial conditions such that the massless fluid flux dominates early on so that 
the situation is effectively that of a single fluid and one has the usual Friedmann-Lemaitre-Robertson- Walker 
(FLRW) spacetime. We find that there is a Bianchi I transition epoch out of which the matter flux dominates. 
The situation is then effectively that of a single fluid and the spacetime evolves towards the FLRW form. Such a 
transition opens up the possiblity of imprinting observable consequences at the specific scale corresponding to 
the transition time. 

PACS numbers: 97.60.Jd,26.20.+c,47.75.+f,95.30.Sf 



I. INTRODUCTION 



The cosmological principle states that the Universe is 
homogeneous and isotropic. Given the increased quality 
of cosmological observations, this fundamental principle 
is now becoming testable, and indeed questionable. That 
questions abound in this area is obvious from the fact 
that we do not have a good handle on the nature of dark 
components that dominate the cosmological "standard 
model" [1|. A large number of alternative models and 
theories have been suggested in the literature, but most 
are not particularly compelling. The treatment of the dif- 
ferent matter components, in particular, is often based on 
the notion of coupled perfect fluids or scalar fields. If 
we are to understand the bigger picture, we need to make 
progress on this aspect, especially if we want to be able 
to consider issues like heat flow f2L^, dissipative mech- 
anisms [5 - 71, Bose-Einstein condensation of dark matter 
Em and possibly many others. 
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We argue that a more natural starting point for this 
endeavor would be the relativistic variational multi-fluid 
approach IfTOl that has (so far) mainly been applied to 
neutron star astrophysics fTTl, and recently to relativistic 
beams and shocks |12J. This approach would seem nat- 
ural since there could have been phases during which the 
Universe would have effectively been anisotropic, with 
different components evolving "independently". For the 
most part, models discussed in the current literature, 
including initially anisotropic geometries, describe the 
matter content in terms of either effectively many com- 
ponent single fluid models |Tj |, or a plain single com- 
ponent 114-161; although isotropisation is expected in 
such situations, as required to end up with a realistic 
(read: in agreement with cuiTently available data) model 
1 17 1, interesting new consequences can however be de- 
rived, e.g. by enhancing an initially vanishingly small 
non-gaussian signal Ill8i . 

As an illustration of the fundamental principles in- 
volved, we develop the formalism for determining the 
cosmological solutions to the Einstein equations for 
a general relativistic, two-fluid model coupling matter 
(non-zero rest mass) and "radiation" (zero rest mass). 
Drawing on the experience from other applications it 
would be straightforward to consider other relevant 
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cases, e.g. involving a dissipative heat flow f4l or super- 
fluid condensates |(8] l9l. However, the chosen example 
is perhaps the most conventional, since the leading-order 
thermodynamics of massless particles has some generic 
features (compare, say, a photon and phonon gas), and 
the same for a massive component when the density be- 
comes (relatively) small. 

Within this context, we will demonstrate how the dis- 
tinct fluid motions lead to anisotropy and the spacetime 
metric taking the form of a Bianchi I solution of the 
Einstein equations. This follows since there is a space- 
Uke privileged vector, associated with the relative flow 
between the two components in the problem. It is im- 
portant to understand that, while this feature is natural 
in the multi-fluid context, it can never arise in the of- 
ten considered multi-constituent single fluid. The multi- 
fluid hypothesis implies that each component (labelled 
by an index x) of the matter and radiation sourcing Ein- 
stein equations follows its own timelike vector m^; the 
relative flow between the various fluids then generates a 
privileged spacelike direction along which the Bianchi I 
solution aligns. However, it is important to recognize that 
it is the fluxes ri^ - n^i^, where Mx is the particle number 
density, that are the fundamental sources. In particular, a 
fluid can be moving quickly with respect to another, yet 
if its density is much smaller its flux can be negligible. 

Such a choice is by no means new [19| and recent 
work in given circumstances have shown, here again, the 
possibility of isotropisation |20|, although behavior very 
diff'erent from the standard cosmological one can also be 
found f2T\ . (A useful review on anisotropic solutions 
and their cosmological use is Ref. [22 1.) For instance, it 
has been suggested |23r ,25J that since Bianchi universes, 
seen as averaged inhomogeneous and anisotropic space- 
times, can have effective strong energy condition violat- 
ing stress-energy tensors, they could be part of a backre- 
action driven acceleration model. 

Yet another reason for studying such cosmological 
models stem, curiously, from the observations! Large 
angle anomalies in the Cosmic Microwave Background 
(CMB) indeed have been observed and discussed for 
quite some time [26 -29| and related with underlying 
Bianchi models lf30ll3Tl . It is not our aim here to decide 
whether or not the data do indeed imply some amount of 
anisotropy, but we shall at least assume that they do not 
rule out the possibility altogether. Note in that respect 
that further, currently ongoing observations of different 
backgrounds will determine, for instance, if the CMB 
dipole is fully originating from mere local Earth motion 
(and should thus be removed altogether from the data) or 
if part of it is cosmological ||32| . 

In order to remain close to the observationally ver- 
ifiable model, we shall concentrate on the example of 
the radiation to matter transition for which, in principle, 
the underlying microphysics ought to be well-known, up 
to the a priori necessarily negligible Dark-Matter to ra- 



diation coupling. We then ask whether it is possible 
to have a cosmological epoch where there is a relative 
flow of radiation with respect to the matter, but out of 
which the expansion becomes isotropic and the relative 
flow dissipates. We will demonstrate that the short an- 
swer to this question is yes, as flux domination of one 
fluid over the other leads to an effectively one-fluid sit- 
uation, thus yielding an effective Friedman-Lemaitre- 
Robertson-WaUcer (FLRW) Universe. In essence, the 
cosmological principle appears to be satisfied on both 
sides of the transition, but the transition itself puts for- 
ward a Bianchi I behavior with a spacelike privileged 
direction. Our goal here is to, first of all, establish this 
possibility and then consider the compatibility of such a 
model with current observational data 13311341 . 

On the technical side, the two-fluid nature of the prob- 
lem introduces several terms that are not present in the 
one-fluid case. We will "skew" the discussion some- 
what by introducing variables that were found useful in 
the stability analysis of two-fluid systems by Samuelsson 
et al. Il35l . In particular, we will take into account the 
fact that two-fluid systems have two speeds of "sound", 
and use causality to constrain parameter values that enter 
through the equation of state. We will also introduce the 
so-called cross- constituent coupling, which occurs when 
the equation of state has terms containing both fluid den- 
sities. It is an equilibrium property and thus is non- dis- 
sipative. While the coupling is not the main focus here, 
it is important for a follow-on analysis |36| where we 
consider so-called two-stream instability. This can oc- 
cur when there is a relative flow between two fluids with 
cross-constituent coupling. If a disturbance is developed 
on top of the relative flow, and the coupling is strong 
enough, it can become unstable if it appears to move, 
say, to the right with respect to one fluid, but to the left 
with respect to the other. In this sense, the work here 
has the additional purpose of building the "background", 
relative-flow configurations. 

The outline of this paper is as follows: In Sec. [ll] 
we construct cosmologies having two Killing symme- 
tries, with the subsequent Einstein tensor components 



presented in Sec. II A Sec. II B contains a brief review of 
the two-fluid formalism and how it applies in the current 
context. We also show how our formalism can be im- 
mediately employed to describe relativistic condensates 
(which reduces to the standard descriptions of terrestrial 
systems, such as superfluid helium four). In the follow- 



ing Sec. Ill we show how an ideal gas in the presence of 
a radiation field leads to a system with cross-constituent 
coupling, and then construct a simpler model containing 
similar characteristics. Sec. |IV] restricts the analysis by 
removing the spatial-dependence in the metric and mat- 
ter (The more general set of equations are required for 
the two-stream instability analysis of 1361 .) This same 



section includes a numerical analysis subsection |IV B| 
and ends with a discussion of the results. We finish with 
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some concluding remarks in Sec.|V]and an appendix con- 
taining more details on how the equations are obtained. 



n. COSMOLOGIES WITH TWO SPACELIKE 
KILLING VECTORS 

We will choose the simplest possible two-fluid model: 
the relative matter flow is in one direction (to be taken 
along the z "axis"), and orthogonal to it will be two, 
mutually orthogonal spacelike Killing vector fields (one 
along the x "axis" and another along the y "axis"). We 
will use as our x and y coordinates the two parameters 
that naturally generate the Killing vector fields and 
J/'^. With this choice we have 

^'^ = (0,1,0,0) , J/'' = (0,0, 1,0). (1) 

It is also the case that 



= g.^X''^' = gn- 



(2) 



Finally, if we let t denote the time coordinate then the 
two symmetries imply the remaining metric components 
are fiinctions of only z and t. 

There is some remaining freedom in the choice of co- 
ordinate system, i.e. it can be shown that the so-called 



synchronous gauge (goo = - 1 and go/ = 0) that reduces 
the metric to 



As' 



-df ^ +g,Ax^ +gyAy^ +g-,Az^+2g,AxAz+2gy-AyAz, 

(3) 

can be utilized. Within this gauge choice there is another 
change of coordinates that can be made, namely t - t, 
X - x(x,z), y - y(y,z), and z - z, that sets the terms g\^ 
and g23 to zero. The final form of the metric is thus 



d/ 



-Af- + AIAx^ + A^.Ay^ + A^Az^, 



(4) 



where the As (N = {x, y, z}) are, as yet unknown, func 
tions of f and z- When the z-dependence is relaxed, the 
spacetime described by Q is of the well-known Bianchi 
I type. Although we focus on this case later in Sec. FV" 
we keep the z-dependence here because a follow-on anal- 
ysis 1361 will need the full z-dependent equations. 



A. The Einstein Tensor 

The non-trivial Einstein Tensor coefficients can be 
straightforwardly computed with the known geometric 
quantities given in the Appendix. Letting a dot " ' " and a 
prime " ' " denote, respectively, d/dt and d/dz, we have 



G', = 
G\ = 



- [HMy + HM, + HyH-) + ^ [/;. + /; + (/, + - /,/, - /,/, - lyl, 

- (//,, + H-) - [hI + HyH, +Hl) + ^ [/; + (/,. - /,) ly] , 

[/, + + (H, - H,) /, + [Hy - //,) /,,] , 



(5) 



where we have introduced the "Hubble"-like functions 

(N = {x,3;,z)) 



Ah 
Ax 

and the "inhomogeneity" functions 



5< 

Ah 



(6) 



We will see below that when the z-dependence is 
dropped, the two-fluid energy-momentum-stress compo- 
nents are such that T * v - T^\., implying for the Einstein 
tensor G '^x = G-'y. 

Clearly, not all these components can be independent 
of each other, for otherwise the overall problem would 
be ill-posed because of too many equations. But recall 
that there is the Bianchi Identity VyG''^ = 0, which for 
(7) the situation here yields two independent components: 
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= G', + dfi', + [h, + Hy + //,) G', + [h + ly + I-) G~, - H,G \, - HyG\. - H,G\, 
= G'. + dfi\ + [H^ + Hy + H,)G', + [l_, + Iy)(G'.,-G\,). (8) 

I 



It is important to note that the second of these vanishes 
identically when there is no z-dependence, because then 
the Einstein tensor component - 0. This means that 
we still need three metric degrees of freedom. 

B. General Relativistic Two-fluid Formalism 

We will use the formaUsm developed by Carter ifTOl 
and various collaborators (see Andersson and Comer 
fi 1 1 for a review and references). The fundamental fluid 
variables consist of two conserved number density four- 
currents, to be denoted n'^. Recall that x is a constituent 
index (for which there is no implied sum when repeated). 

From the currents, we can form three scalars, namely 

A so-called "master" function -A(n^,«y,«xy) (the two- 
fluid analog of the equation of state) is assumed, which 
plays the role of Lagrangian for the system. The energy- 
momentum-stress tensor is 

T%^^>6^, + n^,^il+ny,, (9) 

where 

"V^A-^fi^-ny^ (10) 
is the generalized pressure and 

is the chemical potential covector It is also the momen- 
tum canonically conjugate to the current n^- 

Formally, the ^'^^ and S" coefficients are obtained 
from A via the partial derivatives 

dnj^y 5«x 

The fact that the momentum fi^ is not simply propor- 
tional to the corresponding number density current n^' is 
a result of entrainment (as it is known in the neutron star 
literature; see, for example, [37J): the motion of one fluid 
induces a momentum in the other fluid, and vice versa. 
Entrainment vanishes if the coefficient is zero. 

Finally, the equations for each fluid consists of a con- 
servation equation 

= 0, (13) 



and an Euler equation 

= 0, (14) 

where the vorticity two-form is defined by 

co^, = 2V[X]. (15) 

the square brackets indicating antisymmetrization of the 
enclosed indices. It is important to understand that the 
condition V^T'^y = is satisfied once the equations of 
motion are satisfied. Contrary to the single-fluid case, 
^ijT^v - does not yield enough information to com- 
pletely determine the two-fluid evolution. 

Note that the above way of writing each Euler equa- 
tion makes manifest its geometric meaning as an integra- 
bility condition for the corresponding vorticity, a point 
that has been much emphasized by Carter 1 10] (see also 
IfTTI ). It also immediately supplies a formalism for su- 
perfluid condensates, since setting yU^ = VyOx (where 
Ox represents the phase of the relevant quantum wave- 
function) guarantees that the fluid vorticity is zero. The 
non-relativistic limit of the fluid equations in this case 
recovers those that are well-known for, say, helium su- 
perfluids. 

The symmetries do much to simplify the fluid equa- 
tions. It is easy to see that the vanishing of the Lie deriva- 
tive of with respect to X'' and J/'^ requires to be a 
function only of t and z- We also assume that «x is or- 
thogonal to and J/'^. The unit four- vectors take the 
form 

< = (<,0,0,Mi) , ^^-y/r^ZMi)^. (16) 

The entrainment parameter becomes 

«xy = "xWy {uWy ' ^z"x«y) ' (17) 

while the momenta reduce to 

//J = (SXMx + -^"Swy) ■ ^^^^ 
Finally, the components of T^^y are 

T', = -H nxM^yU? + nywX' 
T\, = P'y = 
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(19) 



where 



>P = A - «x {u'ji'^ + Kijf) - tiy [uy, + u;i4) . (20) 

The remaining item required to completely specify the 
matter is a particular form for the master function A. This 



we will provide in Section III 



We see from the above that our problem has been re- 
duced to finding solutions for the four matter variables 
{rix, Mj} and the three metric functions An. The conserva- 



tion equations ( 1 3 1 now take the form 



= ^ [A^AyA,n^u[) + — (AAyA.n^ui) , (21) 
while the Euler equations reduce to 



dz dt 



(22) 



The remaining equations are those of Einstein, con- 
structed from Eqs. (|5]l and ( 19 1. 

As mentioned in the introduction, we introduce some 
new variables that are convenient for the multi-fluid 
analysis. Since there are two fluids we have the well- 
established result of two modes of "sound" propagation 
IflOM : namely. 



(Jlnju" 
Slnwx 



(23) 



These are "bare" in the sense that they only equal the lo- 
cal wave speed when there are no interactions between 
the fluids |35|. A measure of the interactions are the 
cross-constituent couplings defined — slightly modified 
from ll35l — as 



5 In /i" 



d In n„ 



C. 



yx, 



'y 

where, if we set the entrainment to zero. 



(24) 



(25) 



The Cxy represent a key channel through which the two 
fluids "see" each other (especially when the entrainment 
is zero) |35, 36|. 

Some final words on this set-up is about our frame of 
reference. We have chosen a frame that is not attached 
to either fluid. One might expect it would be easier to 
work in either of the fluid rest-frames, but this is actually 
not the case. Starting with the metric in Eq. Q, we can 
show that "jumping" on a fluid rest-frame introduces a 
shift vector into the metric. 



Let xf^ be the rest-frame coordinates of, say, the x- 
fluid. We can assume that the coordinate transforma- 
tion does not involve the orthogonal pair {x,y], so that 
t - t(t, z), X - x,y - y, and z - z{t, z), which guarantees 



Mx - 0. What we want is = 0, which implies 



dz 

dt 



dziq 
dz ui ' 



(26) 



and therefore z must depend on both t and z- We can now 
assume that f = f. However, the change of coordinates 
also affects the metric; in particular. 



dz 

'dt 



(27) 



III. 



A COSMOLOGICAL TWO-FLUID SCENARIO: 
MATTER AND RADIATION 



When the particle species of a fluid has mass m", it 
can be useful to separate out from A mass density terms; 
namely, 



A - -m^n^ - nv'ny - &(n\, n^, ) , 



(28) 



where & contains other information about the fluid ther- 
modynamics, and relative motion effects. The two-fluid 
cosmology we have in mind has a combination of "mat- 
ter", with mass m" = m, and "radiation", which means 
= 0. We assume a non-zero cross-constituent cou- 
pling and zero entrainment. One of our conserved cur- 
rents is the total particle flux of the matter. Since we are 
ignoring dissipation in the flows, we can use the total en- 
tropy flux of the system as our other conserved current. 
The bulk of this is due to the radiation. To simplify the 
notation, we set Wx - n, Uy - s, f/- = jj, and fi^ = T, 
which is the temperature. 

To see how a cross-constituent term can come about, 
consider the usual way of combining a (non-relativistic) 
gas and radiation in the energy density and pressure: 



mn H — riT + aT , 
2 



riT + -ar 
3 



(29) 



(30) 



where a is constant. We take as our fundamental ther- 
modynamic variables n and s, and so the temperature, 
obtained as T = dp Ids, is a function of both. Hence, the 
ideal gas contribution will generate a cross-constituent 



coupling (cf. Eq. (24i). Even if we take the temperature 
as fundamental, there would still be its coupling with n. 

Writing T in terms of {n, s] explicitly is not tractable. 
So for the purpose at hand, it is perhaps clearer to con- 
sider a simpler, algebraic construction where the depen- 
dence is explicit. With that in mind, we will use a master 
function of the form 



A — —m*n — K^s' 



,4/3 



(31) 
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where we have placed a polytropic coupling to the en- 
tropy in an effective mass m* for the matter; namely. 



(32) 



where o-„ > 1, o"s > 1, and t^s are constants. The re- 
maining fluid variables are 



1 



^' = -K.s'''' + (o-„ + o-,-l) {m* - m) n, (33) 
T 



m + (Ta (m* - m) , 

,1/3 



KsS ' + (Ts (m* — m) ■ 



Cnfi — 



0"s 



(34) 

(35) 
(36) 



It is also worthwhile to consider the other direction of 
the evolution, which is that back to the past, where the 
universe contracts and heats up to the point where the 
temperature scale is much higher than that of the mass 
scale. 



IV. HOMOGENEOUS BACKGROUND 



Assuming that the background is only time-dependent. 



then the two matter equations (21 1 and (22i imply (for 
X = {n,s)) 



where 



jUCn = 0-n(0"n - 1) (m* - m) , 

4 



KsS^^^ + CTs (cTj - 1) (m* — m) ■ 



(37) 



(38) 



There are a few comments to be made about this con- 
struction. In order to have a model that cools as it ex- 
pands, we see that n — > and 5 — > which also means 
m* — > m. This also ensures that the "dust" limit of the 
standard cosmological scenario /i m and Cn — > is 
achieved. Finally, we recover the usual result for the 
massless fluid of 5 oc and — > 1/3. In fact, if we 
eliminate the second term in p8|) using ([35]l, we find 



s - crT^ and cr - 



3 (fT, - 1 - cl) 



(39) 



A^/z^V, = , A,AA«xVl + ^x =^x, (40) 



where Alx and AT^ are constants and we have introduced 
Vx = A-u^. One can also show that T~, = is automat- 



ically guaranteed by (40 1, provided that the integration 
constants satisfy 



M„A^n + MsM=0. (41) 
Equation ( [40| allows, in principle, to write «x and Vx in 
terms of Ajs,-, which can be put into the Einstein equations 
to get a closed system of equations. But, since we will 
be solving the equations numerically, it is actually easier 
to use the original differential equations, which can be 
shown to take the form 



1 + y 2 1 + y 2 
. 1 -H y 2 1-1- y 2 > 

for the densities and 



Vn ' 






{ h \ 










( C ^ 
•-n '-'ns 




n 














S 







yvj ^ * ^ 

for the velocities. 

To solve for the metric we use the definition of and three of the Einstein equations (setting = M^j^ y^ith Mp, 
the Planck mass) to evolve {Aj>j,//n) as follows: 

= + H,H, - ^ (1 + 2y2) + Ts (1 + 2y2)l , (44) 





1 


n 




s 




V s ) 





+ H, + 
H, + H, + 



\ + vl 

Hz 

1 4- V2 i 



(42) 
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Hy = -Hl,+HMz-^[t^n{\+2V^^ + Ts{\+2Vl)\, (45) 

PI 

H, = -Hl + H,H,- ^iixn + Ts), (46) 

PI 



The so-called Hamiltonian constraint G', = — ^rT', is 

M2 



+ //,//, + HyH, = ^ (-A + ^„y2 + r,9y2) . (47) 



M2 

PI 

r 



The initial conditions therefore consist of four mat- 
ter, and six metric initial conditions; i.e. the set 
{n(fo), s(fo), V„(fo), Vs(fo),A!s(fo),//!s(fo)), where to is the 
initial time. 



A. Preliminaries: The Matter Quadratures 



It is useful at this point to apply the results of Eq. (40 1. 
For the matter we can write 



1 



n — 



Vn 



A-ji' 



^\+Mll(A,jif 



(48) 
(49) 



The first indicates that the metric coefficients must grow 
with time if both n — > and yu — > m (the conditions for 
cooling); that is, we can be sure that our model allows for 
both expansion and cooling. The second relation there- 
fore shows that Va — » with time. 

A similar analysis for the entropy fluid is complicated 
by the fact that it is massless, and thus the associated 
chemical potential (i.e. the temperature) can go to zero. 
In particular, it is not clear a priori that the entropy fluid 
velocity 



(50) 



remains finite, i.e. whether or not A,T grows with time. 
Actually, the entropy relation 



1 



s - 



(51) 



shows that 5 — » even if A-T — > 0. In fact, we see that 



< A ,AjA,i < N,. 



As for the behavior of A^T, we can show that substituting 
Eq. (HT} into Eq. ([39]) results in a cubic for {AJIM^f, 
which is 



IA,T 

Mm 



A^ 



A Ay I \a-Ml 



(53) 



If the FRLW solution is obtained in the late time limit, 
then the last term tends to a constant, and hence A,T as 
well. This shows that Vj — > const and m~ — > 0. In fact, 
we see in Fig.[T]that this is precisely the case. The bottom 
line is that this form of model is such that the expansion 
can become isotropic, and damp out the three-velocities 
of each fluid. 



B. Numerical Results 



Numerically solving the system of Eqs. (42 1-(46 1 re- 
quires that we rewrite those in terms of dimensionless 
quantities. Rescaling the time variable to f — > m^f/Mp,, 
and denoting by an overdot the derivative with respect to 
this new dimensionless time, we set 



« 



X = 



m-' 

together with 
yielding 



y = 



ft = — and 
m 



T = 



T 
m' 



T = Tns'W 



3(iT„-o-s)-4 



P = 1 + cr„fx"»-Vr, 



(54) 

(55) 

(56) 
(57) 



showing that the system is fully determined provided the 
two arbitrary dimensionless constants a-j and f are given. 

To make comparison with standard cosmology clearer, 
we further rewrite the Bianchi I metric Eq. (HI in the form 



(52) di^ = -dfHa2(f)(e^^'dxHe2^'d/ H-e^^'-dz^), (58) 
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FIG. 1: Velocities derived from Eq. l |43[ >, as functions of the 
e-fold number N (defined in the main text). The underlying 
parameter values for this plot are a-j = 1, f = 0.1, iXn = 1.1, and 
as = 1.1. The initial values are such that «(0) = 3.9 x 10"', 
s(0) = 1, V„(0) = 0.99, VM = -4.25 x 10"^ As(0) = 2, 
and HxiO) = 2.89 (for each N). This figure illustrates that Kn, 
although initially very large, rapidly decays to zero while Vs 
remains essentially negligible, and almost constant, at all times. 



with ^ /3s - 0, thus defining the scale factor a(t). The 



relations to pass from Eq. Q to Eq. (58 i are then 



A^AyA., and /5x 



In 



3 AyA^ 



(59) 



with similar relations for yS, and obtained by circular 
permutations of the indices (x,y,z)- The so-called shear 
variables ||T| are given by 



Al 



(60) 



We can rewrite the equations of motion in terms of the 
e-fold number A^, defined through 



a(t) = e'^, 

by using the relation 

d 1 / v d 

— = - [Hj, + Hy + HA — 



dA^ 



(61) 



(62) 



The figures that illustrate our results all use this parame- 
ter for the horizontal axis. 

A realistic model having two FLRW phases connected 
by a Bianchi I transition is realized through numerical 
solutions of Eqs. (|42| -(46 1. We use the exact solutions 
of Eqs. ([48] l and (|40|l, together with the Hamiltonian con- 
straint ( |47| l as a measure of the numerical eiTor. This is 
given in Fig. |2j which shows the relative error, for our 
particular choice of parameters, to be limited to at most 
10-1^ 



J 




• Hamiltonian constraint 

■ n — n^x^i 



FIG. 2: Evaluation of the overall numerical errors associated 
with the solution of Eqs. |42]l-|46]l: with the same parameters 
as in the previous figures, we show here the levels at which 
Eqs. l |48[ l, \5i) and \A1) are satisfied. The most error-prone sit- 
uation, at the end of the calculation, still satisfies the constraints 
to better than 10"'^ 

Fig. |3] shows the behavior of the fluid variables with 
around the radiation to matter transition, i.e. with the 
state parameter w smoothly varying from its initial value 
of J to zero. The rescaled number density x is found to 
be negligible throughout, even though its contribution to 
the energy density eventually dominates. The rescaled 
entropy y provides, roughly, all of the energy density p 
initially and for most of the transition, but eventually be- 
comes negligible, as expected. Finally, the temperature 
T is seen to decay to zero, while the rescaled chemical 
potential /} asymptotically takes its fiducial value unity. 

The behavior of the metric and the shears are displayed 
in Figs. |4] and |5] respectively. The beta coefficients 
change from being initially equal to each other (zero in 
the numerical calculation), to final constant values. With 
a rescaling of the spatial coordinates we can absorb these 
constants so as to return to the usual FLRW metric. Here 
we use the same parameters as before, except that we 
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■ A (f = 1) 

13, (f = 1) 

- A (f = 1) 

- A = 10) 

■ A P = 10) 

- A (f = 10) 



3x10" 

o 

^2x10-*^ 
IxlO-f* 
0x10-'' 




FIG. 3: Background fields for the same parameters as Fig. 
[T| The thick solid line corresponds to the state parameter w, 
i.e. the ratio of p (thin full line) and p = (thick double-dotted 
line), the number density .i; oc n is essentially negligible at these 
scales, while the density p is dominated by the contribution of 
the entropy y s (thick dashed line); the conjugate variables, 
namely the rescaled temperature T (thick dotted line) and mat- 
ter chemical potential /i (thick long dashed line) both decrease, 
with /( — ^ 1 as expected. 



have taken f = 1, 10. The reason is illustrated in Fig.|5] 
which shows that the shears cr^, initially vanishing (be- 
cause we start with a FLRW radiation dominated phase), 
increase first during the transition, reach a maximum and 
eventually decrease to vanishingly small values, which 
is expected for the final FLRW matter dominated epoch. 
As one might expect, as the coupling f is increased, the 
anisotropics increase. 

Finally, Fig. [6] shows the time evolution of the sound 
speeds and c^, as well as Cns and Csn- We find that 
c\ - 1/3 throughout, although c\ is modified in the radi- 
ation era. Both cross-coupling terms decrease with time. 
This, along with the vanishing of the relative flow, is an 
important result for a two-stream instability analysis, for 
it implies that the conditions for instability are naturally 
eliminated by the overall expansion of the universe. 

With the solution at hand, it is now possible to return 
to the original equations and understand what is taking 
place during the transition. Originally, we set initial con- 
ditions in the radiation era, for which T s » pm, with 
the extra requirement that Vs <*; 1 : this means that the 
evolution of the three Hubble functions, and hence of 
the scale factors, will be identical, so the shears vanish 
and we are in a FLRW phase. Then, as the product [in 
begins to grow, with T s decreasing, the matter velocity. 



FIG. 4: Metric coefficients j8, as functions of the e-fold number 
N, for the same parameters as Fig. [T| except that f = 1, 10. 
The full thick line represents /?;, while the other two (dotted 
and dashed lines) stand for ji^ and yS,. respectively, satisfying 
Px = Py = - in agreement with our setting. 

provided it was large enough to begin with (and we see 
numerically that we need to set it very close to unity in 
order to have a visible effect) is still large enough that the 
corresponding term becomes important and the scale fac- 
tors begin to evolve in different ways. Finally, even this 
velocity becomes sufficiently small with respect to unity 
that one recovers the FLRW symmetry expected for the 
matter dominated epoch. As to why the matter velocity 
can be large, and yet the universe be radiation dominated, 
is because it is the flux that enters the Einstein equations; 
a large velocity can be compensated for by a small num- 
ber density. 

V. CLOSING REMARKS 

A main goal of this work was to develop for cosmol- 
ogy the general relativistic, multi-fluid model (derived 
from a variational formalism) that has so far been used 
mostly for neutron star astrophysics. While we consid- 
ered only a two-fluid system, the formalism itself can, 
in principle, handle a number of different fluids. As it 
comes from an action, coupling to other fields can be im- 
posed in more or less standard ways. For example, elec- 
tromagnetism can be incorporated via the usual gauge 
coupling, thus allowing for plasmas and their effects on 
the system. We also demonstrated how relativistic con- 
densates follow automatically because the formalism is 
written in terms of the conjugate momenta, and simply 
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FIG. 5: Shear functions cr, as functions of the e-fold number 
A', for the same parameters as Fig. [T] except that f = 1, 10. 
The full thick line represents cr-, while the other two (dotted 
and dashed lines) stand for cr,. and cr,, respectively. Because of 
the relation between the metric coefficients it turns out that 
the shears satisfy a similar relation, namely cr^ = cTy = -jcr-. 
We also see that the anistropies increase as the coupling f is 
increased. 



FIG. 6: The bare sound speeds c^, cj, and cross-constituent 
couplings C„s and Cs„ as functions of the e-fold number A', for 
the same parameters as Fig.[T] Since = 1/3 throughout, it is 
not shown in the figure. Both cl (dotted line) and C„s (dashed 
line) are smoothly decaying functions of time, while Cm (full 
line, hardly visible on the figure) is essentially negligible at all 
times. 



setting the momenta to be gradients of scalars automati- 
cally insures zero vorticity. 

The two-fluid model we introduced is valid for ap- 
plications in cosmology. Even though the relative mo- 
tions were anti-aligned, the model illustrated behavior 
that one might expect from a close examination of the 
radiation-to-matter transition. The main task was to build 
a model in which both the radiation and the matter dom- 
inated eras were describable by means of an FLRW met- 
ric, as necessary to fit nucleosynthesis, CMB, and large 
scale structure formation data ll33l l34l . We found that 
such a situation could easily be implemented, provided 
the relative fluid velocity is large enough at the transi- 
tion time, an assumption that needs to be justified on the 
basis of primordial cosmology models. Indeed, in the 
now well-established framework of inflation B8l , it is 
tremendously difficult for the Universe to remain with 
any relevant amount of left-over anisotropy: in fact, in- 
flation was precisely invented to, among its expected out- 
sets, remove any primordial anisotropy! 

It should be recalled at this point that the model pre- 
sented above is the simplest set-up of what might be en- 
visioned for the transition itself, for we have not taken 
into account, for example, the non-conservation of the 
photon number through its coupling with luminous mat- 
ter, matter flows with more than one constituent, or rela- 



tive flows at arbitrary angles. Obviously, one would not 
be too surprised if comparison to observational data indi- 
cated the need for a more elaborate model. Note also that 
we have assumed here the matter fluid to have only one 
flux-component. This may not be a reasonable assump- 
tion; it is, however, largely a scale-dependent statement. 

There are not so many ways to produce such primor- 
dial anisotropy. Among the most natural are perhaps 
models based on some amount of non trivial electromag- 
netic phenomena taking place during very early epochs. 
Consistent with large scale astrophysical observations of 
y-ray halos around active galactic nuclei |39|, the ex- 
istence of relatively intense intergalactic magnetic fields 
of the order of 10 G have been deduced, whose for- 
mation is expected to be of primordial origin. Some 
inflationary models POl are able to produce such large 
scale magnetic fields, that are statistically isotropic. It re- 
quires a special effort to construct a so-called "hairy" uni- 
verse Bn in which the resulting magnetic field (or any 
other gauge field coherent over large distances) points 
to a privileged spatial direction; off-diagonal TB and EB 
spectra could be induced by such models B2l . hence pro- 
viding an observational means to validate them. 

A special spatial direction can also exist in more rad- 
ical scenarios. In one such model, for instance, a pla- 
nar domain wall remains all through the inflation phase, 
thereby breaking the rotational invariance of the final 



11 



perturbation power spectrum 1*431 . Multifield inflation 
can also serve that purpose by producing vorticity, al- 
though at second order in the relevant perturbations Il44l . 

Another way to induce a non FLRW universe (per- 
haps the simplest) is to start with a theory having a built- 
in privileged timelike vector with which the dominant 
fluid may not necessarily align; examples are provided 
by the Hofava-Lifshitz setup [45, 46 1, originally aimed 
at renormalizing gravity, and the Einstein-aether theory 
§1-7 1.' Depending on the initial conditions, their solution 
can actually also relax to the usual FLRW solution, the 
fluid unit vector then aligning dynamically with the aether 
vector 149115011 . 

In all these situations, it remains to be seen whether 
some cosmological variables might take values that diff'er 
from their canonical ones, as derived in the framework 
of the best-fit vanilla single field inflation paradigm. To 
clarify the situation, a full perturbation theory should 
now be examined |51|. Unlike those that assume only 
a single fluid, a perturbation analysis of a two-fluid sys- 
tem has to take into account the possibility of two-stream 
instability |35|. 

Such instabilities are well-established in plasma 
physics, and have also been argued for in laboratory 
superfluids and their neutron star analogs l52l l53l . 
Samuelsson et aJ. |35 1 have shown that a relative velocity 
and some type of coupling (cross-constituent or entrain- 
ment) for a system of two relativistic fluids is a necessary 
condition for two-stream instability. Roughly, there is a 
"window" of instability that opens when a mode appears 
to be, say, right-moving with respect to one of the flu- 
ids, yet left-moving with respect to the other In this 
paper we have shown that the conditions for such in- 
stabilities to exist, a relative flow between two coupled 
fluid components, may be satisfied in cosmology. We 
have also shown that cosmological expansion provides 
a mechanism for shutting down the instability by clos- 
ing the window, since both the relative velocity and the 
cross-constituent coupling are driven to zero. 

If such instabilities were to be triggered, a basis for a 
set of observational constraints (or possible detections) 
for the transition epoch may be established. In a com- 
panion paper (36"], we explore whether these instabili- 
ties develop before the instability window is closed. If 



an instability were to develop in some of the anisotropic 
transitions, they would most definitely leave relevant im- 
prints in both CMB and large scale structure data, in the 
form either of non-gaussianities, bizarre polarization dis- 
tributions and spectra, and special scales corresponding 
to the Hubble volume at the transition time. For instance, 
it can be argued that such instabilities can occur during 
the matter to cosmological constant transition if and only 
if the latter is made of a fluid, hence having a state pa- 
rameter w > -1; however close w is to -1, such a fluid 
could initiate an instability that an actual cosmological 
constant, having w - -I, could not. Therefore, observ- 
ing the relevant consequences of these instabilities at the 
relevant length scales would allow a discrimination be- 
tween these two otherwise indistinguishable models. 
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Appendix: Geometric quantities 

For the metric given in Q we find the Christoff^el co- 
eflicients to be 



r' 

.v.v 








^ tx 




r?,, = //v , r 




xz 


= Ix , 






r 

XX 


-(^: 







leading to the following non-vanishing components of 
the Ricci tensor: 



//, J]h^ + //; - ^ [/; + Ix [ix + ly - Iz)] , 
//, ^ Hk + h; - -1 [/; + /, (/, + /, - /,)! , 
= H,Y^H^ + Hi-^ [/; + /; + /, (/, - /,) + /, (/, - /,)] , 



(64) 
(65) 
(66) 
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and scalar 



(67) 
(68) 

(69) 



with the sign convention for the Riemann tensor given by not all independent. 

From these, one can obtain the Einstein tensor (|5]l. As 
pointed out in the main text, the Einstein equations are 



[1] P. Peter and J.-P. Uzan, Primordial cosmology (Oxford 
Graduate Texts, Oxford University press, UK, 2009). 

[2] B. Modak, J. Astrop. Astr 5, 317 (1984). 

[3] J. Triginer and D. Pavon, Class. Quant. Grav. 12, 689 
(1995). 

[4] N. Andersson and C. S. Lopez-Monsalvo, Classical and 

Quantum Gravity 28, 195023 (201 1). 
[5] S. Weinberg, Astrophys. J. 168, 175 (1971). 
[6] L. K. Patel and S. S. Koppar, Australian Mathematical 

Society Journal Series B - Applied Mathematics 33, 77 

(1991). 

[7] H. Velten and D. J. Schwarz, JCAP 1109, 016 (201 1). 
[8] P. Sikivie and Q. Yang, Physical Review Letters 103, 

111301 (2009). 
[9] T. Harko, Phys. Rev. D 83, 123515 (2011). 
[10] B. Carter, in Relativistic Fluid Dynamics (Noto, 1987), 
edited by A. Anile and M. Choquet-Bruhat (Springer- 
Verlag, Heidelberg, Germany, 1989), vol. 1385 of Lect. 
Notes Math., pp. 1-64. 
[11] N. Andersson and G. L. Comer, Living Reviews in Rela- 
tivity 10, 1 (2007). 
[12] E. Nakar, A. Bret, and M. Milosavljevic, Astrophys. J. 
738, 93 (2011). 

[13] A. Gromov, Y. Baryshev, and P. Teerikorpi, Astron. As- 
trophys. 415, 813 (2004). 

[14] A. Emir Giimrukfiioglu, C. R. Contaldi, and M. Peloso, 
JCAP 11, 5 (2007). 

[15] C. Pitrou, T. S. Pereira, and J.-P Uzan, JCAP 4, 4 (2008). 

[16] H.-C. Kim and M. Minamitsuji, Phys. Rev. D 81,083517 

(2010) . 

[17] P. -P. Dechant, A. N. Lasenby, and M. P. Hobson, Phys. 

Rev. D 79, 043524 (2009). 
[18] A. Dey and S. Paban, ArXiv e-prints (201 1). 
[19] P. Sandin, General Relativity and Gravitation 41, 2707 

(2009). 

[20] T. Harko and R S. N. Lobo, Phys. Rev. D 83, 124051 

(2011) . 

[21] S. Calogero and J. M. Heinzle, Physica 240, 636 (2011). 



[22] C. G. Tsagas, A. Challinor, and R. Maartens, Phys. Rep. 
465, 61 (2008). 

[23] J. D. Barrow and C. G. Tsagas, Classical and Quantum 

Gravity 24, 1023 (2007). 
[24] K. S. Adhav, S. M. Borikar, M. S. Desale, and R. B. Raut, 

EJTP 8, 319(2011). 
[25] M. Cataldo, E Arevalo, and P. Mella, Astrophys. Space 

Sci. 333, 287 (2011). 
[26] D. J. Schwarz, G. D. Starkman, D. Huterer, and C. J. Copi, 

Physical Review Letters 93, 221301 (2004). 
[27] C. J. Copi, D. Huterer, D. J. Schwarz, and G. D. Starkman, 

Advances in Astronomy 2010, 847541 (2010). 
[28] L. Perivolaropoulos, ArXiv e-prints (201 1). 
[29] Y.-Z. Ma, G. Efstathiou, and A. Challinor, Phys. Rev. D 

83, 083005 (2011). 
[30] A. Pontzen and A. Challinor, MNRAS 380, 1387 (2007). 
[31] A. Pontzen, Phys. Rev. D 79, 103518 (2009). 
[32] D. J. Eixsen and A. Kashlinsky, Astrophys. J. 734, 61 

(2011). 

[33] E. Komatsu et al., Astrophys. J. Suppl. 192, 18 (2011). 

[34] W. J. Percival, B. A. Reid, D. J. Eisenstein, N. A. Bah- 
call, T. Budavari, J. A. Erieman, M. Eukugita, J. E. Gunn, 
Z. Ivezic, G. R. Knapp, et al., MNRAS 401, 2148 (2010). 

[35] L. Samuelsson, C. S. Lopez-Monsalvo, N. Andersson, and 
G. L. Comer, General Relativity and Gravitation 42, 413 

(2010) . 

[36] G. L. Comer, P. Peter, and N. Andersson, to be submitted 

(2011) . 

[37] G. L. Comer and R. Joynt, Phys. Rev. D 68, 023002 
(2003). 

[38] M. Lemoine, J. Martin, and P. Peter, eds.. Inflationary 

Cosmology, vol. 738 (2008). 
[39] S. Ando and A. Kusenko, Ap. J. Lett. 722, L39 (2010). 
[40] M. M. Anber and L. Sorbo, JCAP 10, 18 (2006). 
[41] M.-A. Watanabe, S. Kanno, and J. Soda, Physical Review 

Letters 102, 191302 (2009). 
[42] M.-A. Watanabe, S. Kanno, and J. Soda, MNRAS 412, 

L83 (2011). 



13 



[43] C.-H. Wang, Y.-H. Wu, and S. D. H. Hsu, ArXiv e-prints 
(2011). 

[44] A. J. Christopherson, K. A. Malik, and D. R. Matravers, 

Phys. Rev. D 79, 123523 (2009). 
[45] R Horava, Phys. Rev. D 79, 084008 (2009). 
[46] D. Bias, O. Pujolas, and S. Sibiryakov, Phys. Rev. Lett. 

104, 181302 (2010). 
[47] T. Jacobson, PoS QG-PH, 020 (2007). 
[48] T. Jacobson, Phys.Rev. D81, 101502 (2010). 
[49] W. Donnelly and T. Jacobson, Phys. Rev. D 82, 064032 



(2010) . 

[50] 1. Carruthers and T. Jacobson, Phys. Rev. D 83, 024034 

(2011) . 

[51] T. G. ZIosnik, ArXiv e-prints (2011). 

[52] N. Andersson, G. L. Comer, and R. Prix, Mon. Not. R. As- 
tro. Soc. 354, 101 (2004). 

[53] N. Andersson, G. L. Comer, and R. Prix, Phys. Rev. Lett. 
90,091101 (2003). 



